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Abstract 

As a continuum work of Bhaumik et al who derived the common eigenvector of the number-difference 
operator Q = (a'a — P&) and pair-annihilation operator ab (J. Phys. A9 (1976) 1507) we search for 
the simultaneous eigenvector of Q and (ab — a'b') by setting up a complex differential equation in the 
bipartite entangled state representation. The differential equation is then solved in terms of the two- 
variable Hermite polynomials and the formal hypergeometric functions. The work is also an addendum 
to Mod. Phys. Lett. A 9 (1994) 1291 by Fan and Klauder, in which the common eigenkets of Q and 
pair creators VP are discussed. 


1 Introduction 

Coherent states [1-2] are widely used in many aspects of quantum physics. The bosonic coherent state, 
which describes the quantum state of laser, obeys coordinate-momentum minimum uncertainty relation and 
possesses non-orthonormal and over-complete properties. Generalized coherent states has been constructed 
and applied by theoreticians since 1970s, among which the coherent state | q,a) for charged bosons [3] is a 
remarkable one, when one introduces ’’charge” by defining two types of quanta possessing ’’charge” +1 and 
-1 with corresponding annihilation operators a and b, so the operator Q = a)a — b'b is endowed with charge 
operator, | q,a) is constructed based on [Q,ab] = 0, [a,V] = [6, 6+] = 1. In the Fock space the charged 
coherent state is 


I q,a) 
Q \q,a) 


cX 


I n + q,n), 


n=o \J(n + q)\n\ 

q\q,a), ab\q,a) =a\q,a) 


where C q is the normalization constant. In quantum optice theory \q. a) was named by Agarwal as pair- 
coherent state in [4], Q is the two-mode photon number-difference operator. By observing [Q, W] = 0, in 
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Ref. [5] Fan and Klauder also constructed the common eigenvector of Q and a'b^ with use of the Dirac’s 
^-function in countor integral form proposed by Heitler [6]-[7]. One then naturally to ask waht is the 
simultaneous eigenstates, denoted as \q, k ), of Q and (ab — This question is full of physical meaning 

in quantum optics, because most nonlinear interactions in tha parametric approximation reduce to a bilinear 
form 

Hi = itm (ab — a^V) ; (1) 

since 

[(a6 - a l '& t ) ,Q] =0, (2) 

where k is related to the susceptibility in the parametric process. Hi is responsible for producing two-mode 
squeezed states [8-11], which is not only useful for optical communication and weak signal detection, but also 
embodies quantum entanglement, i.e. the Einstein-Podolsky-Rosen (EPR) correlations [12] for quadrature 
phase amplitude are intrinsic to two-mode squeezed light, the idler-mode and the signal-mode generated 
from a parametric amplifier are entangled each other in a frequency domain. Solving this question is quite 
difficult. We recall Dirac’s guidence [13]: “When one has a particular problem to work out in quantum 
mechanics, one can minimize the labour by using a represenation in which the representatives of the more 
important abstract quantities occurring in that problem are as simple as possible”. At first glance, we 
thought that the charged coherent state representation | q, a) was a good candidates for tackling the problem 
as simple as possible, but after some tries we found that it was not. Eventually we find that the entangled 
state representation [14-15] is of assistence for searching for the desired common eigenvector | q, k) of Q and 
[ab — cR&i), 

Q\q,k) = q\q,k) , (3) 

(ab - a f 6 f ) | q, k) = {q - k - 1) | q, k) . (4) 

because it possesses well-behaved features. Thus we shall briefly review the properties of entangled state 
|£) in Sec. 2. Then in Sec. 3 we shall make full use of |£) to set up a complex differential equation for the 
overlap (£| q , k). In Sec. 4 we employ the two-variable Hermite polynomials and the hypergeometric function 
to solve the differential equation. The Gauss’ contiguous relation of hypergeometric function is essential for 
us to derive the common eigenvector of Q and (ab — . 


2 Brief review of the bipartite entangled state representations 

In [14] and [15] we have constructed the bipartite entangled state 

|0) = |£), £=|0e*, (5) 


exp 


id 2 


-■^-Ka t +O t -a t & t 


|£) satisfy the eigenvator equations 

(a + b')\o = t; 10, (a f + fe) I0=f 10- 


( 6 ) 


\ (Xi + X 2 ) 10 = 10, (Pi - PO 10 = ^2 10 , (7) 

i.e. |0 is just the simultaneous eigenvector of two-particles’ center-of-mss posititon X\ — X 2 and the relative 
momentum Pi — P 2 in Fock space, we name it the EPR eigenstate since EPR were the first who used the 
commutative property [Xi — X 2 ,P\ + P 2 ] =0 to challege the incompleteness of quantum mechanics [11]. 
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According to Dirac’s representation theory [13]: “To set up a representation in a general way, we take a 
complete set of bra vectors, i.e. a set such that any bra can be expressed linearly in terms of them,” we use 
the normal ordering form of the two-mode vacuum state projector 

| 00 ) < 00 | =: e - a]a ~ bH :, ( 8 ) 

and the technique of integration within an ordered product (IWOP) of operators [16]-[17] , we can prove 




3 Complex differential equation for the new state \q,k) 

By noticing (5) we see 

a 10 = (C-^)IO, 6|0 = (r-a f )l0, 

(0<* f = (0(0-0, (06t = (^|K -a), 

(0 (ab - a f 6 f ) = (C| [ab - (C* - 6) 6 f ] = (C| [&C - (£ - a) 0 + 1] • 

«*io = (!4)k>, ‘’ie-(^ + i)ie. 


so 


Then using 

we re-write (11) as 


<£|(a&-at&+) = (0 




ac 2 

Operating (13) on the state \q,k), and using (4) we obtain a complex differential equation 


(Cl (ab~a)b ] ) \q,k) 

" d * d 

C^+C ^r + 1 


(Cl Q,k) = (q-k- 1) (Cl q,k) 


On the other hand, from (5) we have 

(a'a-b'b) |C) = (Ca + - & + C*) 10 

= ICI ( e~ lv> a 1 — e* v &l) exp 

- 


10 s 


+ |C|(e-^a t +e iv & t )-a t 6 t 


100} 


which together with (3) leads to another differential equation 

(Cl (a 1 a - b ] b) \q, k) = (C| q, k } = q\q, k) 

In the next section we want to solve (14) and (16). 


(9) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


(15) 


(16) 
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4 Solution to Eq. (16) 


After many tries we find the solution of (14) and (16) is 

(£1 <1 , k) = £(k)e~^ 2/2 A, (17) 

where C(fc) is the normalization constant, which keeps undisturbed in the following calculations, so we neglect 
it in the derivation process, and 


OO 7 

A = J 2 (£*•0 2*1 (- n , 2 + 1;9 + 1 ; 2 ) 

n—0 


2 .Pi is the hypergeometric function defined as 


ir■ ( a /3- 7 . e ) = V' ( Q )» (^)n "" = r (^) y' r(n + a)r(n + / 3 ) 2 


the symbol (a) means 


„=o Wn n! r ( a ) r (/ 3 )^ ) F(n + 7) n!’ 


(«)„ = F = a (a + 1) (a + 2) • • • (a + n - 1), 


r(a) 

and the two-variable Hermite polynomials is defined as [18] 

min/m -kO 


Hm,n&0= E 


mini 


l=o 


l\{m — l)\(n — l)\ 


(-1 ) l c-T n -= Hn ,m(V,0, 


whose generating function is 


°° 


E — T—rHm,n(€, C) = exp (-ft ; + ft) + f'£*) • 
^' mini 


m,n=0 


(18) 


(19) 

( 20 ) 


( 21 ) 


( 22 ) 


Note that the convergent condition for the hypergeometric function defined in (19) is \z\ < 1, 7 7 ^ 0, —1, —2, •• 
•, so 2 -Pi(— n, | + l;g+l; 2 )isa formal power series. 

Now we prove (14): Firstly, we notice 


(§- ( (SI «, k) = S + e- KI ’ /2 |/l) , (23) 

and 

f ^ (£1 Q, k) = r (-| e_l€|2/2 ^ + e- |5|2/2 ^^) , (24) 

it then follows 

+ r (Cl <7, fc) = -ICI 2 <£| q , fc> + e-l«l 2 / 2 + r ^ (25) 

so 

(^ +r ^ + 1 ) (cl9,fc) (26) 

= e_l?|2/2 (^+r^ + l-|^| 2 )E^W^>^ 2p 1 1 (-n,^ + l; g + l;2). 
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(27) 


Then we use the property 


^ H m ,n r,o = nH m)n _ x (r,o, (c,o = mH m _i n (r,o 


+ nH m 

,n —1 — £ ^m,n) -^m,n+l + m,H m -1 

,n — 

which can be derived from (21) and (23), as well as 

|£| 2 -ffm >n = c (H m+ltn + nH ra , n _i) 

= 33m+l,ro+l + U777 H rn —I—1 + (777 + 71+1) -£7yn,n> 

we have 

(^Jr + r^ + i) (eu,fc) 

= e _l «| 2 / 2 X)^Lo ^[C^n+g^-l +C(<? + «)^n+«-l,n 
33ro-t-g-|-l,rt-t-l 77/(71 + q)H n -\-q—\, n — X (g + 2?l) 3Tn+(7,n] 2-^1 ( 77, ^ + 1, g + 1, 2) 

= e“l«| 2 / 2 X)(T=0 ^ {« [#n+9,n + («■ + g)-ff„+ g -l,„-l] 

+ (g + 7l) [H n +q jn + TlH n ^.q — \ n —\\ 

l,n+l 77.(77 + Q^H n -\-q— l,n—1 (g + 2 ti) 2 -Ci( 77, ^ + 1, g + 1, 2) 

= e“^l 2 / 2 J2^=o ^i[(q + n ) nH n+ q -i,n-i - H n+ q + i, n +i} 2^1 (-n, § + 1; g + 1; 2) 
= e _ ^l / 2 X)^ 0 ^,i7„+ g ,n[(g + 77 + I) 2 ^i (—n — 1, | + 1; g + 1; 2) 

—7i 2 -Fi( — ^ + 1, f + 1; g + 1; 2)]. 

Then using two Gauss’ contiguous relations about the hypergeometric function [19-20] 

2 F 1 (a,f3-,r,e) = 2^1 (P, a; 7 ; e), 


[7 - 2/3 + (/3 - a)e] 2 Fi (a, /3; 7 ; e) + /?(1 - e) 2 ^i (a, /3 + 1; 7 ; e) 
-(y - P) 2 F 1 (a,/3 - 1; 7 ; e) =0, 


and letting a = | + 1, /3 = — 77 , 7 = g + 1, e = 2, we have 

{q-k - 1) 2 Cl (— + 1, — 77 ; g + 1; 2) + 77 2 -Fi(— + 1, —77 + 1; g + 1; 2) 
— (g + 77 +I) 2 Ci(— + 1, —?7 — 1; g + 1; 2) = 0 
Substituting (30) into the right-hand side of (30) we finally obtain (recovering €(k)) 

{^ +C w +1 ) (el q,k) 

00 1 h 

= (q-k- l)£(k)e-M 2/2 J2 -,Hn +q ,u (f, O 2 M-n, - + 1; g + 1; 2) 

n=0 

= (q-k-l)(£\q,k), 

thus we have proved the solution of (14). On the other hand, from (21) we know 

H n+q , n (£*,£) = e~ iqv H n+qtn (|£|, |£|), 

so the solution (14) automatically satisfied with (16). The solution seems new. 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
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5 The simultaneous eigenstate of Q and (ab — a)b t) 

Now we hope to obtain \q, k) from the information of (£| q, k ). Using the completeness relation (9) of |£) 
and (17)-(18) we can have 

\<l,k) = J ^ I 0 ( 0<?,0 (36) 

/ _j2 £ 050 1 ht 

— 10 e-'« |2/2 E ~ H n +q ,n (0, 0 2 fi(-n, - + 1; q + 1; 2). 

77 n! 2 

Then using (22) we expand |£) in (5), 

10 = e -'« |2/2 £ ^^(00)10,0) (37) 

1,3 =0 ' J ' 

OO 

= e -'« |2/2 E -m= 7 H ij(£,C)\l,j), 

where |Z,j) is the two-mode Fock state. Substituting (37) into (36) and using the integration formula 

J ^e-l e\ a H lt j ( 00 ) = S l!m 6j, n m\n\, (38) 

we have 

°° 1 7 

|g, fc) = £(fc)E-r o +1;9+1;2 ) ( 39 ) 

' n! 2 2 

n=0 

/ t2 <7 00 -1 

—e" 1 ^ E - m Hu(U*)HZ +g ,n(te*)\U) 

= 00 E — E_ g ^' 2 0(~n, ^ + 1; q + 1; 2 ) |n + g, n) 

n—0 V n! 2 

050 7"j"Tl L, 

= £(fc)a^E^4r 2Fl(_n ’f + 1;9 + 1;2)|0i0) - 

n—0 

As we have known that the convergent condition for the hypergeometric function defined in (19) is \z\ < 1, 
7 /O, —1, —2, • • •, so \q 1 k) seems not normalized as a finite number, but as a divergent one. To see this more 
clearly, using (17)-(18), (9) and (38), formally we have 

(q,k\q,k) = J ^ (q,k\ 0 (£\q,k) (40) 

/ j2 £ 1 J 

— E 1 H n +q ,n (0,0 2^1 (-n, - + 1; q + 1; 2) 

77 ^0 n! 2 

00 

x E 

n'=0 

= |C(fc)|-XI hL±-2l! | 2 F,(-n i+ l ; , + l;2)| 2 . 

77-. Z 

n=0 


2jtfn'+ g ,n' (0 0) 2 -Fi(-0, | + 1; q + 1; 2) 
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Therefore, the common eigenvector of Q and ( ab— , like the common eigenvector of Q and aftf, is 

normalized as a singular function, so its applications are greatly limited. However, the exploration for formal 
solution of |< 7 , k) has its own sense in mathematical physics. 

In summary, as a continuum work of Bhaumik et al [3] we have set a complex differential equation in the 
entangled state representation for deriving the new the common eigenstate \q, k) of parametric interaction 
Hamiltonian and number-difference operator, the complex differential equation has been solved in terms of 
hypergeometric function and the two-variable Hermite polynomials and the solution seems new. Thus this 
paper embodies a new use of hypergeometric functions and is of theoretical mathematical physics meaning. 
This work is also an addendum to Ref. [5], in which the common eigenkets of <5 and two-mode creators a'b' 
are discussed. 
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